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Abstract 

The perturbation of the dromion of the Davey-Stewartson-1 equation is studied over 
the large time. 

1 Introduction 

In this work we construct the asymptotic solution of the perturbed Davey-Stewartson-1 
equations (DS-1): 

1 "32 , »2> 



id t Q + + di)Q + (Gi + G 2 )Q = eiF, 

(1.1) 

%Gi = ~d v \Q\ 2 , d v G 2 = -^dAQ\ 2 . 

Here e is small positive parameter, a = ±1 correspond to so called focusing or defocusing 
DS-1 equations. 

The equations (1.1) describe the interaction between long and short waves on the liquid 
surface in case the capillary effects and potential flow are taken into account [||, gj. The 
existence theorems for the solutions of this equations in the various functional classes are 
well-known jl . The inverse scattering transform method for the DS equations was 
formulated in |J-[f|. This method allows to construct solitons |J and to study global 
properties of the generic solutions, for instance, the asymptotic behaviour over large times 
[10, [ll]]. The asymptotics for nonintegrable DS equations were studied in [12]. 



The perturbations of the equations ( |1.1[ ) arise due to a small irregularity of bottom or 
by taking into account the next corrections in more realistic models for the liquid surface 
considered in 0, §|. For the first case the perturbation takes the form: F = AQ. Here A 
is real constant and its sign corresponds to decreasing or increasing depth with respect to 
spatial variable £. 
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We start with the solution of DS-1 constructed in ||: 

p(erit . p) = pA/xexp(it(A 2 + /x 2 )) 

2cosh(/xe)cosh(Ar ? )(l-a|p| 2 ^(l + tanh(Ar ? ))(l+tanh(/ie)))' 

where A, p are positive constants defined by boundary conditions as r\ — > — oo and £ — > 
—oo; p is free complex parameter. 

This solution decreases with respect to spatial variables exponentially. Besides if a = 1 
and ^-\p\ 2 > 1 this solution has singularities at some lines r] = const and £ = const. 

The solution ( |1.2|) was called dromion in work Q. The inverse scattering method for 
dromion-like solution of (1.1) was developed in Jjjj. Prom the results of @, and |L0| it 
follows that the soliton of the DS-2 equations is unstable with respect to small perturbation 
of the initial data. These results stimulate the studies of the dromion perturbation. 

To construct the asymptotic solution over the large times we suppose the parameter p 
depending on a slow time variable r = et and obtain an explicit formula for the p(r). 

The possibility of the full investigation of the linearized DS-1 equations plays the main 
role to construct the perturbation theory of the nonlinear DS-1 equations. For (1+1)- 
models it has been found in [|nj, that spectral functions of the Lax pair form the basis of 
solutions of the linearized nonlinear PDEs. The same idea is true for the DS-1 and DS-2 
equations, as shown in [16], [17| and it is used below. 

The asymptotic analysis given in this paper is valid for the solutions ( |1.2|) without the 



singularities. It means if a = 1, then ^-\p(t)\' 2 < 1. If the coefficient of the perturbation 
A > 0, then \p\ increases with respect to slow time. It allows to say that the singularity 
may appear in the leading term of the asymptotics as r — > oo. However we can't say this 
strongly for perturbed dromion in our situation, because our asymptotics is usable only 
when r <C log(log(e -1 )). In general case the appearance of the singularities in the solution 
of nonintegrable cases of the Davey-Stewartson equations is known phenomenon [18|. 

The obtained result for the problem about the interaction between the long and short 
waves on the liquid surface shows when the depth decreases the formal asymptotic solution 
can be described by adiabatic perturbation theory of the dromion for the focusing DS-1 
equations at least for |t| <C s^ 1 log(log(e -1 )). 

The contents of the various sections are as follows. Section 2 contains a statement of a 
problem and a result. Section 3 is a brief treatment of solving of linearized DS-1 equations 
using the basis formed by solutions of the Dirac equation. In Section 4 we apply the results 
of Section 3 and construct the first correction of an asymptotic expansion for the solution 
of the problem formulated in Section 1. Section 5 is devoted to reducing of a modulation 
equation for the dromion parameter. In appendix we demonstrate explicit formulas for 
functions which are used to obtain the basis set solving of the linearized DS-1 equations. 



2 Problem and result 



We construct the asymptotic solution of equations (1.1) on mod(0(£ 2 )) uniformly over 
large t. The perturbation operator is F = AQ and the boundary conditions for G\ and 
G2 are: 



G 



l|£->-oo = U\ = 



2 cosh 2 (At?) 



G 



2|?7- 



U 2 = 



2 cosh 2 (pO 



(2.1) 
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We seek the asymptotic solution as a sum of two first terms of asymptotic expansions: 

77, t, e) = W{£, v, t, t) + eU (£, v, t, r), 

Gi(e,J7,*,e) = gi(t,V,t,T) +eVi(£,ri,t t T), (2.2) 

G 2 (C, v, t, e) = g 2 (C, V, *, r) + eV 2 (£, 77, t, r), 
where r = et is slow time. The leading term of the asymptotics has the form: 

W(£,ri,t,T) =q(Z,V,t;p(r)), 
and gi, g 2 are: 



ded v \Wtf,ri,t,T)\ 2 , 



Denote by 7(7") = 1 — o"^|p(r)| 2 and 70 = 7(0). The final result is given by 
Theorem 1. If 

Arg{p{r)) = const, 



7(r) 



exp(2Ar) 
To ) 



where 70 > 1 at a = —1 and < 70 < 1 at a = 1, then the asymptotic solution (2.i ) with 
respect to mod(0(e 2 )) is useful uniformly over t = 0(e _1 ). 

Remark 1. When the time is larger than e _1 , namely, t <C log(log(e -1 )), the formulas 
( |2.2| ) are asymptotic solution of ( |1.1| ) with respect to mod(o(l)) only. 



Remark 2. D. Pelinovsky notes to author, that the modulation of the parameter \p\ may 
be obtained out of the "energetic equality" [fill : 



dt 



d£drj\Q\' 



d£dr](QF - QF). 



3 Solution of linearized equations 

In the present section we obtain the formulas for solution of the linearized DS-1 equations 
on the dromion as a background: 

id t U + (flf + 3%)U + (Gi + G 2 )U + (Vx + V 2 )Q = iF 

dzV 1 = ~d v {QU + QU), d v V 2 = ^d^QU + QU). ( ' > '' ) 



The results of inverse scattering transform J?] and the set of the basic functions [16] 
are used for solving of the linearized equations by Fourier method. It is reminded in 
Subsection 3.1. 



However in contrast to the work [16| here the solution of the DS-1 equations with 
nonzero boundary conditions is considered. It leads to changing of the dependency of 
the scattering data with respect to time (see also 0) and of the formulas which define 
the dependency of Fourier coefficients of the solution for the linearized DS-1 equations in 



contrast to obtained in [16]. This is explained in Subsection 3.2. 
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3.1 Basic set for solving of the linearized DS-1 equation 

In the inverse scattering transform one use the matrix solution of the Dirac system to 
solve the DS-1 equations (see (5|-0): 




^ z k (3.2) 



l % 

^ o d v 

Let ip + and ip~ be the matrix solutions of the Goursat problem for the Dirac system 
with the boundary conditions (see |5J): 

V>nl^-oo = exp(ikr]), ^k-»-t» = 0, 

^il^— >°o = 0> ^22 1 *7 — >— °° = exp(— ifc£); 

V^nk^-oo = exp(i/cr/), V^k^oo = 0, 

Denote by Wfj) j = 1) 2, the columns of the matrix ip + . The column ip~k-. is the solution 
of two systems of equations. An additional system of time evolution is 



<w+, = + i( J _°j ) o { - a„) 2 «S , 



One can obtain the equation like this for the other columns of the matrices if) . These 
equations are differ from (3.4) by the sign of ik 2 in first term of the right hand side only. 

Below we write two bilinear forms defining analogues of the direct and inverse Fourier 
transforms. First bilinear form is 

OO POO 



(X) /")/=/ / d£,dr){x\Hx of + X2M2 /)• (3.5) 



—00 j —00 



Here Xi anci Mi are the elements of the columns x an d M- 

Denote by 4>{i) and <^ 2 ) the solutions conjugated to rpf^ and tp^ with respect to the 
bilinear form ( |3.5| ). 

Using the formulas for the scattering data (Pj) one can write these data as: 

Sa (fc,0 = ^Wfofov^Ewiiltyq, (3.6) 

s a(M) = ^{'>P(2)(€,rj,k),E(2)(ilri))Q- (3-7) 

Here -E'(z) = diag(exp(z), exp(— z)). 

It is shown (0) the elements of the matrices tp^ are analytic functions with respect to 
the variable k when ±Im(k) > 0. Using the scattering data one can write the nonlocal 
Riemann-Hilbert problem for the ipT^ and ipf 2 011 the real axes ([@]): 

^h(C) 7 ?)^) = exp(ifcry) + exp(i/c77)( exp(-ikrj) / d/ si(fc, l)ipf 2 (£,, rj, I) 



^tifaVity = exp(-ifc^) ( exp(i^) / 

V J — ! 



dls 2 (k,l)^ u (^,rj,l) 
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Here 

± 



_ i r */(*') 



2ivr 7_ 00 k' - (k ± iO) 
The Riemann-Hilbert problem for and ^22 nas the form: 



^2l(^V,k) = exp(i/cr?)l exp(-i&r?) / dl si(k,l)i>J 2 (£,r],l) 

V J —00 

^22(^V, k ) = exp(-«'fcf) +exp(-i/c^)( exp(i/c£) / dl s 2 (k, 1)^(^,7], I) 
Introduce second bilinear form 

/*00 /'OO 

(*,//>.= / / dkdl(x\l)^(k)s 2 (k,l) + X 2 (l)» 2 (k) Sl (k,l)), (3.8) 



00 j —00 



where x J is the element of the row x- 

Denote by <p^\ j = 1,2, the row conjugated to = ( <0+] with respect to the 



bilinear form ( |3.8| ). Formulate the result about the decomposition obtained in [16]. 

Theorem 2. Zei Q be such that d a Q G L\ n C 1 /or |a| < 3, if a function f is such that 
d a f G ii n C 1 /or |a| < 4, then one can represent the f in the form 

f = zl (ij,P)(Z,n,l),<pW(Z,V,k))f, (3-9) 

7T 7 



w/iere 

/=^(^5)(c.^fc)^(i)(e,»/,0)/- (3-io) 

3.2 Solving the linearized equation using Fourier method 



In the preceding subsection we have state that the bilinear forms ( |3.5| ) and (3.8) may be 
used to decomposition like as a Fourier integrals. The present subsection shows how to 
use this decomposition for solving of linearized DS-1 equation. 

Theorem 3. Let Q be the solution of the DS-1 equations with the boundary conditions 
—00 — ui and G2I77— ►— 00 

= U2, and Q satisfies the conditions of Theorem 2, the 
solution of the first of the linearized DS-1 equation is smooth and integrable function U 
with respect to £ and rj, where d a U G Li n C 1 and d a F G L\ n C l , for \a\ < 4 and 
t G [0,T ]. Then 



8 t U = i(k 2 + l 2 )U + / dk'U(k-k',l,t)x(k')+ / dl'U(k, I — I', t)n{l') + F. 
J— 00 J — 00 

(3.11) 

If the boundary conditions for the solution of the DS-1 equation equal to zero, then 
X = k = 0. In this case the formulas of the Theorem 3 allow to solve the linearized DS-1 
equation in the explicit form. It was be done in [16]. In contrast of [16], we consider 
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here the solution of the DS-1 equation with nonzero boundary conditions. It leads to 
integral terms in the formula ( 3.11| ). In order to solve the linearized DS-1 equation we 



must transform the formula ( 3.1 1| ) . In the right hand side of ( p. 11 ) the integral terms 



are the convolutions. Go over to the equations for the Fourier transform of the functions 
U(k,l,t,r) with respect to variables k and I. As a result we obtain the linear Schrodinger 
equation: 

idtU + (of + 3 2 )U + M£m) + Mv^W = F. (3.12) 

Here 

U(£,rj,t,r) = — / dkdlU(k,l,t,r) ex.p(—ikr] — 

If , ( 3 - 13 ) 

F(£,r),t,T) = — / dkdlF(k,l,t,r) exp(—ikr) — 
27r Jr 2 

The same equation without the right hand side (as F = 0) was obtained in j7|] for the time 
evolution of the scattering data for the DS-1 equation. 

One can construct the solution of the Cauchy problem for the equations (|3.12j ) with the 
initial condition U(£,rj) = by the separation of the variables. In our case the solution of 
the equations ( |3.12| ) obtained by the Fourier method has the form: 

U(£,rj, t) = — / dmdnU(m, n, t)X(£, m)Y(rj, n) exp(— it(m 2 + n 2 )) 
2n Jr 2 

+ I dnU ll (n,t)Y(r ] ,n)X fl ^)exp(-it(n 2 - n 2 )) 

+ -^= I dmU x {m,t)X(£,m)Y x (r])exp(-it(rn 2 - \ 2 )) 



2vr 

+ U^xX^OYxiv) exp(it(^ 2 + A 2 )). (3.14) 
Here we use notations: 

pttanh(^) +im , . . 1 
X(m,£) = : exp(-im£), XJQ- 



im — [i ^ ^ 2cosh(^£) ' 

Y(n,rj) = - + exp(— inrj), Y\(rj) — 



in — A ' 2cosh(A?y) ' 

d t U = i(m 2 + n 2 )U + F(m, n, t), dtU^ = i(n 2 - /j, 2 )U + F M (n, t), 

d t Ux = i(m 2 - \ 2 )U + F x (m,t), 8 t U^ = -i{\ 2 + fi 2 )U + F^t), (3.15) 

U\t=o = U^\t=o = U\\t=o = Ufj,\\t=o = 0; 



F(m,n,t)= / d£dnF(£,ri,t)X(m,£)Y(n,Tj), 

JR 2 

F ll (n,t)= [ dtdvF&n^X^OYin^), 

JR 2 

F x (m,t,r)= f d£driF(Z,Ti,t)X(m,Z)Y x (Ti), 
Jr 2 

F„x(t,r)= [ d^dvFiCri^X^OYxiv)- 
Jr 2 



(3.16) 
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To solve the linearizing equations we must apply the formulas ( 3.10| ), ( |3.13D , ( [3.16 ), 
solve the Cauchy problems ( |3.15| ) and apply formulas ( [3.14 ), inverse Fourier transform to 
U and then apply formula ( |3.9| ) for U. As a result we obtain the solution of the linearized 
DS-1 equations. The functions Vi,2 may be obtained by direct integrating of second and 
third equations from (|3.1|). 



4 Equation for the first correction 

This and next section contain a proof of the Theorem 1. In this part the equation for 
the slow modulation of the parameter p(r) is obtained. This equation is necessary and 
sufficient condition for the uniform boundedness of the first correction of the expansion 
Q over t = 0(e- 1 ). 

Some complication appears when one use the Fourier method from preceding section 
to solve the linearized DS-1 equations on dromion as a background with slow varying 
parameter p(r). In this case the basic functions used for solving the linearized equations 
depend on slow variable r also and they are the basic set of asymptotic solutions with 
respect to mod(0(e)) only. 

Substitute the formula (|2.2| ) into the equations ( |1.1[ ). Equate the coefficients with the 
same power of e. The equations as e° are realized since W,g±,g 2 are the asymptotic 
solution of nonperturbed DS-1 equations. For the first correction we obtain the linearized 
DS-1 equations: 

id t U + (of + d 2 )U + ( 9l + g 2 )U + {V x + V 2 )W = iH 

a - - - a - - (4.1) 

d(;Vi = -~d v (WU + WU), d v V 2 = —d^WU + WU), 

where 

H = AW — d T W. 

Before to use the formulas from the preceding section we reduce the form of the right 
hand side in first of the equations (|4.1[) . In the leading term the parameter p depends 
on the t only. The other parameters of the dromion (A and p) depend only on the 
boundary conditions and have not changes under perturbation of the equations. For more 
convenience we represent p(r) = r(r) exp(ia(r)) where r(r) = \p(r)\ and a(t) = Argp(r). 
The derivation of W with respect to slow variable r can be written as: 

d T W = d r Wr' + d a Wa'. 

Here the derivatives r' and a' are unknown yet. 

Compute the function H. Using the Theorem 2 and formulas for the functions ^+ and 
4> (see Appendix) we obtain: 

H(k, I, t, t) = exp(-it(A 2 + p 2 )) (p(k, I; p) - R(k, I; p)d T p^j , 

where 



P(k, I; p) = exp(it(A 2 + p 2 ))AW, R(k, I; p) = exp(it(A 2 + p 2 ))drW . 
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In these formulas we have write the dependence on time in explicit form. It allows to 
remove the secular terms (increasing with respect to t) from the asymptotic solution fl2.2|) 
using modulation of the parameters r(r) and a(r). 
The differential equations for U have the forms: 

d t U = i(m 2 + n 2 )U + exp(-it(A 2 + p 2 ))(P(m, n; p) - R(m, n; p)), 
d t U„ = i(n 2 - p 2 )U + exp(-it(X 2 + p 2 ))(P^(n; p) - R^(n; p)), 
d t U x = i(m 2 - X 2 )U + exp(-it(A 2 + p 2 ))(Px(m; p) - R x (m; p)), 
dtU^x = -i(\ 2 + P 2 )U + exp(-it(A 2 + p 2 )){P^x{p) - RM). 

Solutions of these equations are 

p 2 + \ 2 ( V, w \ 

U (m, n, t) = — ^ — ^ -r I P(m, n; p) - R(m, n; p) I exp(-it(p 2 + A 2 )), 



n 2 + m 2 + p 2 + A 2 

U\(m,t) = —5- — j Px(m;p) - R\(m;p) exp(-it(p 2 + A 2 )), 
vnr + p z \ J 

a 2 + A 2 ( » \ 
U^(n,t) = re2 ~ X 2 {P^( n iP) ~ R^( n ;p)j exp(-^0 2 + A 2 )), 

U»\{t) = (P^xim-p) - R^ x {m;p)^t^ V {-it{p 2 + A 2 )). 

One can see that the secular terms may appear because of the last term in the equation 
for Ufj,x- The eliminating of this term leads us to the equation for p{r): 

R/j,X ~ P„X = 0, p| r =o = pa- (4.2) 

As a result U^ x = 0. The another solutions U(m,n,t), U^(n,t) and U\(m,t) are 
bounded with respect to all arguments and over all times. 

One must return to the original of the images U to say about boundedness of the 
solution U,V\,V2 for the equations ( [4.11) . One can see the direct (from U into U) and 
inverse (from U into U) integral transforms as the Fourier transform from the smooth 
and exponentially decreasing functions with respect to the corresponding variables. The 
Fourier transform moves such functions into analytic functions near the real axis. The 
inverse transform moves these analytic functions into the exponential decreasing functions. 
So the solution of ( fO| ) is bounded and decreasing exponentially with respect to the spatial 
variables. 



5 Modulation equation 

Here the equation ( [4.2j ) for the parameter p(r) is reduced to the more convient form. Write 
the derivative of the leading term with respect to the slow variable r. 

r' 2iW r' 

id T W = -a'W - iW- + r . 

r l-or 2 ^(l + tanh(/i^))(l + tanh(Ar ? )) r 
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Denote by T = ^-r 2 and compute the images (-) m a °f ever y term. 

(WU = ^°p(-^' + "')) (ff r - i) log |i - ,r|. 

Denote the image of the last term by h^x. Its image has the form: 

V = ~ g (1 " cT) l yT -^ - 1 - log |1 - aT\ 

The image (•) of AW has the similar form. 

Substitute these formulas into ( |4.2| ) and separate the real and the imaginary parts of 
this equation, then 

ol = 0, 

-(oT-l)logll-oTl --( oT-fl-aniogll-oT 
r r \ 

+ A(aT - 1) log |1 -aT\ = 0. 
Use the notation for T, then the second equation has the form: 

^- = -2<jA(1 - aT) log II - aT\. 

CLT 

This equation defines the evolution of the absolute value of the complex parameter p. The 
argument of this parameter do not change under the perturbation F = AQ. 
Denote by 7 = 1 — err, rewrite the equation for 7. Then we obtain: 

7 ' = 24 7 log(7). 
The solution for this equation has the form: 

7(7-) = exp(Cexp(24r)), 
where 7| r =o = ex P(C)) then we can write the 7(1") in the form: 

/ \ exp(2Ar) 

7(r)=7 • 
The Theorem 1 is proved. 
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Appendix 

A.l Explicit formulas 

Here we remain the explicit forms of the solution for the Dirac equation with the dromion- 
like potential. These forms were obtained in In our computations we use first column 
of the matrix ip + only. 



N roo dp\exp(ikp)) 

iPl! \ _ ( eX.p{lKT]) \ Jri 2cosh(Ap) 



+ 



ni J V / (l-a|p| 2 ^(l + tanh(^))(l+tanh(A7?))) 

/ -o-|p| 2 A^(l+tanh(^)) \ 
x - 8 T h ^L 2^ • (A.l) 

\ 2cosh(^) / 

A. 2 Solutions of conjugated equations 

Here we write the problems for the functions conjugated to ijj^ with respect to the bilinear 
forms. 

The matrix </> is the solution of the boundary problem conjugated to the solutions of 



the problem (|3.2| ), ( f2.1|) with respect to the bilinear form (3J3). First column of the matrix 
4> is the solution of the integral equation: 



1 f v 

^ J-oo 

-1 f 00 



The explicit formula for first column of the matrix <p used in Section 4 for the dromion 
potential has the form: 

f°° dpfiexp(ilp)) 
hi \ _ f exp(i/£) \ J£ 2 cosh(/ip) 



+ 



21 7 V / (l-a| / o| 2 Tf(l+tanh(^))(l+tanh(Ar ? ))) 

( -o-|p| 2 A^(l+tanh(A>7)) \ 
- ap\ exp" A 2 +M 2 ) ) ) • ( A - 2 ) 
2cosh(Ar?) / 



Second bilinear form ( |3.8| ) allows to write the integral equations conjugated to the 
integral equations which were obtained from the nonlocal Riemann-Hilbert equation for 
the matrices ip^ in 0]. Using the equation for the matrix <p(£,T],l) one can show that ipn 
and (f2i are the solutions of the boundary problem for the Dirac system or the equivalent 
integral equations: 

1 p 

<Pll(£,V,l) = exp(?Z£) + - / dr]'Q(C,rj',t)^2i(C,v'J), 

^ J-oo 

^ J -oo 

The functions 9911 and </?2i have the similar form as the functions cpn and 021- 
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